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Abstract. We describe TRILEGAL, a new populations synthesis code for simulating the stellar photometry of any 
Galaxy field. The code attempts to improve upon several technical aspects of star count models, by: dealing with 
very complete input libraries of evolutionary tracks; using a stellar spectral library to simulate the photometry 
in virtually any broad-band system; being very versatile allowing easy changes in the input libraries and in the 
description of all of its ingredients - like the star formation rate, age-metallicity relation, initial mass function, and 
geometry of Galaxy components. In a previous paper (Groenewegen et al. 2002, Paper I), the code was first applied 
to describe the very deep star counts of the CDFS stellar catalogue. Here, we briefiy describe its initial calibration 
using ElS-deep and DMS star counts, which, as we show, are adequate samples to probe both the halo and the 
disc components of largest scale heights (oldest ages). We then present the changes in the calibration that were 
necessary to cope with some improvements in the model input data, and the use of more extensive photometry 
datasets: Now the code is shown to successfully simulate also the relatively shallower 2MASS catalogue, which 
probes mostly the disc at intermediate ages, and the immediate solar neighbourhood as sampled by Hipparcos 
- in particular its absolute magnitude versus colour diagram -, which contains a somewhat larger fraction of 
younger stars than deeper surveys. Remarkably, the same model calibration can reproduce well the star counts 
in all the above-mentioned data sets, that span from the very deep magnitudes of CDFS (16 < i? < 23) to the 
very shallow ones of Hipparcos {V < 8). Significant deviations (above 50 percent in number counts) are found 
just for fields close to the Galactic Center (since no bulge component was included) and Plane, and for a single 
set of South Galactic Pole data. The TRILEGAL code is ready to use for the variety of wide-angle surveys in the 
optical/infrared that will become available in the coming years. 



1. Introduction 

The number counts of Galactic stars in a given bin of ap- 
parent magnitude [mx,m\ + dm a] - where A stands for 
a passband - and towards an element of galactic coordi- 
nates {i, b) and solid angle dfl, is given by the fundamental 
equation of stellar statistics (see Bahcall 1986 for a review) 

/>oo 

N{mx,i,b) ^dmx dr p{r) (j){Mx,r) (1) 
Ja 

where r is the line-of-sight distance, and p(r) is the stel- 
lar density as a function of the position r — {i,b,r). r, 
when measured in parsecs, is related to the absolute and 
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apparent magnitudes Mq^a and mx, and to the interstellar 
absorption Ax, by 

Mo, A = TOA - 5 logr - Ax{r) + 5 . 

(/)(Mo,A,r) is the intrinsic distribution of stellar absolute 
magnitudes, i.e. the intrinsic luminosity function (LF) of 
the stars considered at r. 

To describe the stellar densities p(r) for the largest 
possible volume, and to a lesser extent also 0(Ma, r), is the 
ultimate task of the so-called Galaxy star count models. 
To achieve these goals, the usual way is to assume the 
functional forms of p and (f>, and then compare the results 
of Eq. to observed number counts in several Galaxy 
fields. A number of assumptions help in simplifying the 
task. The first one is to recognize that the Galaxy can be 
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separated in a few distinct components, such as the disc, 
halo, and bulge: 

P = Pd + Ph + Pb , (2) 

each one of these components having a simple expression 
for their density. The second one is to assume an intrinsic 
LF (f) which is virtually independent of r, i.e. (f){M\,r) = 
<p{M\), for each component. 

The recipes for 4>{M) can be of two types. Either one 
(1) assumes an empirical (f){M), derived from e.g. star 
counts in globular clusters or in the Solar Neighbourhood, 
or (2) assumes a theoretical 4'{M), derived from a set of 
evolutionary tracks together with suitable distributions of 
stellar masses, ages, and metallicities. 

Option (1) was the preferred one in the past, and the 
one adopted in some of the most successful Galaxy mod- 
els. Despite their success in reproducing several sorts of 
data, it is not difficult to find points of inconsistency in 
many of such models. A common approximation, for in- 
stance, has been that disc stars of different absolute mag- 
nitude present difi'erent scale heights (and hence Pd(r)): 
Bahcall & Soneira (1980, 1984; and many authors later 
on) assigned the scale height of 325 pc for My > 5.1, of 
90 pc for My < 2.3 dwarfs, and linearly interpolated in 
between. This separation was interpreted as a coarse sepa- 
ration into "young" and "old" populations. Red giants in- 
stead were assigned 250 pc. Gilmore & Reid (1983) adopt 
similar approximations for main sequence stars. Mendez 
& van Altcna (1996, 1998) do the same, and moreover as- 
sume a unique scale height for all evolved stars (subgiants, 
giants, and white dwarfs). From the point of view of stellar 
and population synthesis theories, these approximations 
are clearly not justified, for a series of reasons: (a) Most 
coeval stellar populations contain both red and blue stars 
characterized by initial mass values which are, at least for 
the most luminous objects, very much the same; it is then 
very unlikely that the spatial distribution of these red and 
blue stars could be different, (b) Similarly, young stellar 
populations contain both bright and faint main sequence 
stars, whose relative scale heights cannot change that dra- 
matically in time-scales of less than one Gyr. (c) Moreover, 
it is remarkable that population synthesis theory indicate 
that in star-forming galaxy components (e.g. in the thin 
disc), red giants arc relatively young - most having less 
than say 2 Gyr (see Girardi & Salaris 2001) - and not old 
objects; applying the largest scale heights to all giants is 
then simply wrong. 

Although this sort of inconsistency is not inherent to 
method (1), they are completely removed by the use of 
method (2). In the latter, at any r the relative numbers 
of stars with different colours and absolute magnitude 
strictly obey the constraints settled by stellar evolution 
and population synthesis theories; hence, p{r) cannot be 
arbitrarily changed as a function of absolute magnitude. 
On the other hand, method (2) allows p{r) to be eas- 
ily expressed as a function of other stellar parameters, 
such as age and metallicity - something not possible with 



method (1) where individual stellar ages and metallici- 
ties are not available - then allowing the simulation of 
important effects like metallicity gradients, scale lengths 
increasing with age, etc. This turns out to be a significant 
advantage of method (2) over (1). 

Models that follow method (2) may be put under the 
generic name of "population synthesis Galaxy star count 
models", and have been developed in the last decades 
by e.g. Robin & Creze (1986), Haywood (1994), Ng et 
ai. (1995), Castellani et al. (2002), and Robin et al. 
(2003). These works benefit from the releases of extended 
databases of stcillar evolutionary tracks to predict the 
properties of stars of given mass, age, and metallicity. 
Some assumptions then are necessary to give the distribu- 
tions of these stellar parameters. Such distributions may 
be derived, for instance, starting from an initial mass func- 
tion (IMF), an age- metallicity relation (AMR), and a law 
for the star formation rate (SFR) as a function of Galaxy 
age. 

In the present paper, we will describe a Galaxy model 
developed according to the population synthesis approach, 
taking particular care in the consistency among the dif- 
ferent sources of input data. It has been developed with 
a primary task in head, which is, essentially: to be ca- 
pable of simulating the same sort of data that will be 
released by some major campaigns of wide-field photome- 
try conducted these years. Of primary importance in this 
context, are the several parts of the ESO Imaging Survey 
(EIS; Renzini & da Costa 1997), the Two Micron All Sky 
Survey (2MASS; Cutri et al. 2003), and the Sloan Digital 
Sky Survey (SDSS; York et al. 2000). Present and future 
data from HST deep fields, VIMOS, VISTA, UKIDSS, 
GAIA, might be considered as well. Moreover, our model 
should also be able to take advantage of the extraordinary 
constraints provided by the astrometric mission Hipparcos 
(Ferryman et al. 1997). Of course, a program which meets 
these aims can be applied to any other sort of wide field 
data as well. 

Before proceeding, let us briefly summarize our pri- 
mary objectives and how these translate into technical 
requirements. 

First of all; our primary goal is to simulate the ex- 
pected star counts in several passband systems, such as 
those used by Hipparcos, EIS, 2MASS, SDSS, etc. For do- 
ing so, wc should be able to consistently predict the stellar 
photometry in a lot of different photometric systems. The 
way out to this problem has been settled in a previous 
work (Girardi et al. 2002), that describes a quite general 
method for performing synthetic photometry and deriving 
bolometric corrections from an extended library of stellar 
spectra. Such tables are now routinely produced for any 
new system that we want to compare the models with. 

The second requirement is of being able to simulate 
both very shallow - but of excellent quality - photomet- 
ric data samples as Hipparcos (Ferryman et al. 1997), 
and very deep ones such as the EIS Deep Public Survey 
(e.g. Paper I). The former case implies that we include 
all important evolutionary sequences, such as most of 
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the main sequence, and giants both in the red giant 
branch (RGB) and red clump, which make the bulk of 
the Hipparcos colour-magnitude diagram (CMD). In the 
latter case, we should also include an extended lower main 
sequence, reaching down to visual absolute magnitudes as 
faint as My ^ 30, which corresponds to stellar masses of 
~ 0.1 Mq- Moreover, old white dwarfs start to become 
frequent at such faint magnitudes as well. Therefore, the 
libraries of stellar data should be extended to the intervals 
of very low masses, and to very old white dwarfs. 

It is also clear that these requirements imply, neces- 
sarily, that we opt for the population synthesis approach. 
In fact, there is little hope that we could collect empirical 
data for such a variety of stars, in the several photomet- 
ric systems involved, and with good enough statistics that 
the intrinsic CMDs could be constructed with reliability. 
Theoretical data, instead, is available for all of our pur- 
poses, as will be shown in what follows. Such theoretical 
data is also routinely submitted to stringent tests against 
photometric data, such as the Hipparcos CMD, star clus- 
ters, eclipsing binaries, red giants with measured diam- 
eters, etc. In general, the errors detected in the models 
amount to less than a few tenths of a magnitude, and just 
for some particular stars and/or passbands. Certainly, the 
time is ripe for completely relying in theoretical data in 
Galaxy star count models. 

The plan of this paper is as follows: Section |21 
and 1^1 detail the code and its input data, respectively. 
Sect. 0] describes its initial calibration as performed for 
Groencwegen et al.'s (2002) work, based mostly on EIS 
and DMS data. Sects . l57!^ and I^T^ details and presents the 
fine-tuning of the initial calibration, including additional 
comparisons with 2MASS and Hipparcos data. Sect. 
draws a few final comments and summarizes the main re- 
sults of the present paper. 

2. The code 

We describe here all the necessary input for computing 
a Galactic Model. Normally, this means simulating the 
photometric properties of stars located towards a given 
direction (£, b). This task is performed by the newly devel- 
oped code TRILEGAL, which stands for TRIdimensional 
model of thE GALaxy^ 

The code is written in C language. Its core is made of a 
few subroutines that efficiently interpolate and search for 
stars of given mass, age, or metallicity, inside a database 
of stellar evolutionary tracks. They deal with all the in- 
trinsic properties of stars - luminosity, effective tempera- 
ture, mass, metallicity, etc. These subroutines, developed 
in Girardi (1997), have so far been used in a series of works 
dealing from the construction of theoretical isochrones 
(e.g. Girardi et al. 2000; Salasnich et al. 2000) to the simu- 
lations of synthetic CMDs for nearby galaxies (e.g. Girardi 
et al. 1998; Girardi 1999, Girardi & Salaris 2001; Marigo 

^ TRILEGAL is also a word commonly used to say "very 
nice" in southern Brazil. 



et al. 2003). Another set of routines, more recently devel- 
oped, deal with all aspects related with synthetic photom- 
etry, i.e. the conversion between intrinsic stellar properties 
and observable magnitudes. They rely on the same simple 
formalism described in Girardi et al. (2002). 

2.1. Scheme 

A general scheme of the code is provided in Fig^ It makes 
use of 4 main elements: a library of theoretical evolution- 
ary tracks, a library of synthetic spectra, some parameters 
of the detection system, and the detailed description of the 
Galaxy components. The libraries of evolutionary tracks 
and spectra can be pre-processed in the form of theoreti- 
cal isochrones and tables of bolometric corrections, so as 
to reduce the number of redundant operations during a 
simulation. These are to be considered as "fixed input", 
but can be easily changed so as to consider alternative sets 
of data. 

The instrumental setup specifies, among others: (1) 
The set of filters-|-detector-|-telescope throughputs in 
which the observations are performed; any change of them 
implies the recalculation of the bolometric correction ta- 
bles; (2) The effective sky area to be simulated; the num- 
ber of simulated stars scales with this quantity. 

The several Galactic components (halo, thin and thick 
disc, bulge, etc) are specified by a their initial distributions 
of stellar ages and metallicities (SFR and AMR), masses 
(IMF), space densities, and interstellar absorption. This is 
done separately for each component. The space densities 
are in the form of simple expressions containing a few 
modifiable parameters, to be specified in Sect. 13 below. 

2.2. The simulation and output 

A run of TRILEGAL is formally a Monte Carlo simulation 
in which stars are generated according to the probability 
distributions already described. Eq. ^ is used for predict- 
ing the number of expected stars in each bin of distance 
modulus. For each simulated star, the SFR, AMR and 
IMF are used to single out the stellar age, metallicity, and 
mass. Finally its absolute photometry is derived via inter- 
polation in the grids of evolutionary tracks (or isochrones) , 
and converted to the apparent magnitudes using the suit- 
able values of bolometric corrections, distance modulus 
and extinction. 

During the simulations, a lot of different stellar param- 
eters can be kept in memory and printed out: initial and 
current mass, age, metallicity, surface chemical composi- 
tion and gravity, luminosity, effective temperature, core 
mass, etc. In the case of thermally-pulsing (TP-) AGB 
stars, this information is also used to simulate the pulse 
cycle variations (see Marigo et al. 2003). 

The calculation initially produces a "perfect photomet- 
ric catalogue" , which perfectly reflect the input probabil- 
ity distributions but for the Poisson noise. This catalogue 
can be later degraded by using the known photometric er- 
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Fig. 1. A general scheme of our 
codes. The continuous arrows refer 
to steps which are performed inside 
the TRILEGAL main code and sub- 
routines; they lead to the simulation 
of perfect (i.e. without errors) pho- 
tometric data. The dashed lines refer 
to some optional steps, usually per- 
formed with external scripts, like e.g. 
those mentioned in Paper I, to gener- 
ate catalogues with errors. 



rors, such as photon noise, saturation and incompleteness 
for a given instrumental setup. This second task does not 
belong to TRILEGAL, but is performed by separated sub- 
routines, like for instance the ones described by Paper L 

3. The input datasets 

There are essentially 5 different input datasets in 
TRILEGAL: 

1. tables of stellar evolutionary tracks, that give their ba- 
sic properties (bolometric magnitude Mboi, effective 
temperature TeS, surface gravity g, core mass, surface 
chemical composition, etc) as a function of initial mass 
Mi, stellar age r, and metallicity Z. 

2. tables of bolometric corrections BC\^ for the several 
filters pass-bands Ai, as a function of Tcs, logg, and 
[M/H], as well as the relative absorption in the several 
passbands with respect to V, Ax J Ay; 

3. the IMF 

4. the star formation rate as a function of age, tp{t), and 
age-metallicity relation, Z{t), for the different Galaxy 
components; 

5. the geometry of the Galaxy components, i.e. the stellar 
density p{r) and differential V-band extinction dAv{r) 
as a function of position r. 

They will be discussed separately in the following. 
3.1. Evolutionary tracks 

Based on our previous work on simulations of synthetic 
CMDs for Local Group galaxies, we have assembled a 
large, quite complete, and as far as possible homogeneous 
- in terms of their input physics - database of stellar 
tracks. They are illustrated in the HR diagram of Fig. |2] 
1) For masses between 0.2 M© and 7 Mq, we use the 
tracks from Girardi et al. (2000) , that range from the zero- 
age main sequence (ZAMS) up to either the end of the 
TP-AGB, or to an age of 25 Gyr for the lowest-mass stars. 
Metalhcities are comprised between Z = 0.0004 ([M/H] = 
-1.7) and 0.03 ([M/H] = +0.2). Another set with Z = 



0.0001 ([M/H] = —2.3) and computed with the same input 
physics (Girardi, unpublished) is included. 

The TP-AGB evolution included in these tracks is 
estimated from a simplified synthetic evolutionary code 
(cf. Girardi & BcrteUi 1998, case of eq. 17 plus eq. 20). 
Although this TP-AGB evolution is very approximated, 
it provides a reasonable initial-final mass relation (see fig- 
ure 2 in Girardi & Bertelli 1998), and hence reasonable 
masses for the white dwarfs to be considered below. The 
maximum mass of WDs attained is about 1.2 Mq for the 
lowest metalhcities, and 0.9 for the highest. On 2002, we 
replaced these simplified TP-AGB tracks for more detailed 
ones computed by Marigo et al. (2003; and in prepara- 
tion), which are based on the new formulation for molec- 
ular opacities by Marigo (2002). 

2) For masses lower than 0.2 Mq, and down to 
0.01 Mq, we include the models for very-low mass stars 
and brown dwarfs with dusty atmospheres from Chabrier 
et al. (2000). This provides us with a main sequence (MS) 
going down to luminosities as faint as log(L/Lo) = —5, 
and to effective temperatures as cold as 916 K. For M < 
0.2 Mq, Chabrier et al. (2000) tracks exist only for solar 
metallicity, which are hence adopted for all metalhcities 
in our models. It is also worth remarking that Chabrier et 
al. tracks evolve on time-scales of some Gyr, and include 
a significant fraction of the pre-MS evolution. 

3) Post-AGB and white dwarfs tracks of a suitable 
mass are attached at the end of the TP-AGB phase of 
all stars of initial mass between 0.6 and 5 Mq. We use the 
post-AGB and PNe nuclei tracks from Vassiliadis & Wood 
(1994) down to logL/L© ~ -1.5. We then shift to the 
WD cooling tracks of Benvenuto & Althaus (1999) with 
CO cores, total mass ranging from 0.5 to 1.1 Mq, a 10~^ 
Mq envelope layer of hydrogen, and non-zero metallicity, 
plus their 1.2 Mq model of zero-metallicity. Unfortunately, 
these tracks end at log(L/LQ) = —5, which implies max- 
imum WD ages of about 10 Gyr. To overcome this prob- 
lem, we have artificially extended the tracks up to ages of 
15 Gyr. Notice that no track is available for WDs between 
1.2 and 1.4 Mq (the Chandrashckhar mass). This is not a 
problem at all, since from our AGB tracks we predict no 
WDs within this range of masses. 
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Fig. 2. HR diagram containing all 
tracks assembled for the solar metal- 
licity. Our database contains similar 
data for 5 other values of metallicity. 
In the electronic version of this pa- 
per, tracks from different sources are 
marked with different colours: Girardi 
et al. (2000; black) for most evolution- 
ary phases of low- and intermediate- 
mass stars, complemented with the 
TP-AGB phase from Marigo et al. 
(2003, and in preparation; magenta), 
massive stars from Bertelli et al. 
(1994; green), very-low mass stars and 
brown dwarfs from Chabrier et al. 
(2000; red), post-AGB and PNe nu- 
clei from Vassihadis & Wood (1994) 
complemented with WD cooling se- 
quences from Benvenuto & Althaus 
(1999; both in blue). 



4) For masses higher than 7 Mq we adopt the same 
tracks as in Bertelli et al. (1994) isochrones, but for Z = 
0.0001 and Z = 0.001 where more recent models (Girardi 
et al. 1996, and 2003, respectively) are used. 

We interpolate among all these tracks in order to derive 
stars with intermediate values of mass, age, and metal- 
licity. All interpolations are performed between points of 
equivalent evolutionary status, as usual in the codes for 
generating isochrones (e.g. Bertelli et al. 1994; Girardi et 
al. 2000). Linear interpolations are adopted, with logm, 
\ogt, and [M/H] being the independent variables. 

The complete set of stellar models for solar metallic- 
ity is plotted in the HR diagram of Fig. |3 Similar grids 
of tracks are available also for metallicities Z = 0.0001, 
0.0004, 0.001, 0.004, 0.008, and 0.03 (limited to M < 7 M© 
in the latter case). 

Finally, we remark that the present stellar database 
corresponds to the "basic set" of isochrones as mentioned 
by Girardi et al. (2002) and available at the web page 
http : //pleiadi .pd. astro . it/ isoc_photsys . 00 .html, 
but for three important improvements: the inclusion of 
post-AGB and white dwarf cooling tracks, the extension 
of very-low mass stars and brown dwarfs down to 0.01 
Mq, and the improved prescriptions for the TP-AGB 
phase. 



3.2. Tables of bolometric corrections and absorption 
coefficients 

Once a star of (L, Toff , [M/H]) is selected by the code, 
its bolometric luminosity Afboi — — 2.51ogL -I- Mboi© is 
converted into absolute magnitudes by means of Mx = 
M\jo\ — BC\. The bolometric corrections BC\ are de- 
rived from a large database of synthetic and empirical 
spectra, according to the synthetic photometry procedure 
throughly descrived in Girardi et al. (2002). Importantly, 
this allows the application to a very wide set of photomet- 
ric systems, provided that we deal with 

— intermediate- to broad-band filter sets (otherwise er- 
rors in synthetic colours become significant); 

— VEGA, AB, or ST magnitudes, or systems in which 
some of the photometric standard stars are also well- 
measured spcctrophotometric standards. 

Regarding the spectral library in use, it is essentially 
the same one described in Girardi et al. (2002), but 
now complemented with white dwarf spectra. Namely, the 
sources of spectra are: 

1. For most stars, the "non-overshooting" set of models 
computed by Castelli et al. (1997) with the Kurucz 
(1993) ATLAS9 code (see also Bessell et al. 1998). This 
set covers Tefi from 50000 down to about 3 500 K (i.e. 
from O to early-M spectral types), and metallicities 
[M/H] from +0.5 to -2.5. 

2. Blackbody spectra for stars with Toff > 50 000 K. 

3. Empirical spectra for M giants (Fluks et al. 1994). 
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4. "BDdustyl999" spectra for all dwarfs cooler than 3 900 
K, and down to 500 K (AUard et al. 2000). 

5. Synthetic DA white dwarf atmospheres from Finley et 
al. (1997) and Homeier et al. (1998), for Tcs between 
100 000 and 5 000 K. 

In addition to the BC values for stars of any 
(Teff, logg, [M/H]), we also compute, for any photomet- 
ric system, the ratio between the absorption coefHcient 
in each filter, and the total absorption in the V band, 
A\/Av- As it is well known (see e.g. Grebel & Roberts 
1995; Girardi et al. 2004), this ratio is not strictly con- 
stant but depends on the spectral type of the star and the 
extinction curve under consideration. 

In the present work we compute the Ax /Ay ratio for a 
G2V star (the Sun) subject to mild absorption {Ay < 0.5 
mag) and following the Cardelli et al. (1989) absorption 
curve with Ry = 3.1. The derived Ax/ Ay quantities^ are 
then applied to stars of all spectral types and reddening 
values, although, formally, they are adequate only for low- 
reddening G2V stars in the case Ry =3.1. This approach 
is adopted just for the sake of simplicity. In alternative, 
it is very easy to implement a more accurate approach to 
the problem, which will be followed in future applications. 

3.3. The initial mass function 

The IMF (j)jn is a crucial ingredient because it determines 
the relative numbers of very-low mass stars, that may 
dominate star counts at visual magnitudes fainter than 
~ 22. We have introduced the IMF in a very flexible way, 
so that it that can be easily changed. In order to be able 
to use star formation rates in units of Afo/yr, our default 
IMF normalization is for a total mass equal to 1, i.e. 



m0,„dm = 1 M(T 



(3) 



Our default IMF is a Chabrier (2001) log-normal func- 



tion, 



cx m exp 



(logm — log mo )^ 



2cr2 



whose parameters are a characteristic mass, mp = 0.1 Mq, 
and dispersion, a = 0.627. 

Other commonly-used IMFs, like segmented power- 
laws (Salpeter 1955; Kroupa 2001) and Larson's (1986) 
exponential form, are also included in the code. 

The IMF as given above refers to the mass distribution 
of single stars. Additionally to them, it is very easy to sim- 
ulate non- interacting binaries in our simulations. When so 
required, we adopt the same prescription as in Barmina 
et al. (2002): for each primary star of mass mi, there is a 
probability /b that it contains a secondary, whose mass m2 
is given by a flat distribution of mass ratios comprised in 



^ Extinction coefficients for the Johnson-Cousins-Glass and 
SDSS photometric systems are tabulated in Girardi et al. 



the interval [by,, 1]. This prescription is particularly useful 
for simulating the binary sequences which are often evi- 
dent in CMDs of open clusters. Typical values of /b and &b 
are 0.3 and 0.7, respectively, which we adopt as a default. 

3.4. Star formation rates and age-metallicity relations 

Each Galaxy component is made of combination of stellar 
populations of varying age and metallicity. In our code, 
their distribution is completely specified by the functions 
SFR, 'ip{t), and AMR, Z{t). Both are given in a single 
input file containing, for each age value, 

— the SFR ip{t) (in units proportional to Mq/jt), 

— the mass fraction of metals Z, 

— the logarithmic dispersion of Z, a{\ogZ). 

Whenever necessary, Z is converted into the logarithmic 
metal and iron contents by means of the approximate re- 
lations 



[M/H] = log(Z/Zo) 
[Fe/H] = logiZ/Ze) 



[a/Fe] 



(4) 
(5) 



(2004). 



where [a/Fe] is the degree of enhancement of a elements 
with respect to scaled-solar compositions, and Zq = 0.019 
as in Girardi et al. (2000). These relations provide [M/H] 
and [Fe/H] values accurate to within ^ 0.03 dex, which is 
good enough for our purposes. 

For this paper, a enhancement is taken into consid- 
eration for the low-metallicity Galaxy components (halo 
and old disc). In these cases, we can safely convert a 
given [Fe/H] into Z by means of Eq. and associate 
to that Z the evolutionary tracks computed with scaled- 
solar compositions using the relations provided by Salaris 
et al. (1993). For metallicities higher than about half solar, 
this approximation is no longer valid and it is preferable to 
use tracks specifically computed for a-enhanced composi- 
tion (see Salaris & Weiss 1998; VandenBerg et al. 2000; 
Salasnich et al. 2000). 

The SFR can only be considered well-known for the 
old Galactic components. In fact, for ages close to 12 Gyr, 
a change of age of 1 Gyr causes just small changes in the 
stellar luminosity function, and hence has a minor impact 
on the simulated star counts. For the disc components, 
the SFR is less constrained. Anyway, the Galactic model 
is also relatively insensitive to the disc SFR, at least as 
long as we are not sampling regions at low Galactic lati- 
tude and/or the Solar vicinity (as will be shown later). In 
general, even more important than the SFR is the choice 
of the AMR Z{t) and its dispersion, that may change the 
position of simulated stars in colour-colour diagrams, and 
cause a significant colour dispersion. 

After these considerations, it is convenient to specify 
what are to be considered "default" SFR and AMR - i.e. 
those used in Paper I, and partially also in this work: 

The SFR was assumed to be constant over the last 11 
Gyr for the disc, and constant between 12 and 13 Gyr for 
the halo. 
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Fig. 3. Intrinsic My vs. B~V diagrams that follow from our default choices of SFR, AMR, IMF, evolutionary tracks 
and spectra, shown for both the disc (left) and halo (right). In both cases, the histograms to the right and top show 
the corresponding luminosity and colour functions, respectively. Each diagram contains about 10^ simulated objects. 



The AMR for the disc was taken from Rocha-Pinto 
et al. (2000). [Fe/H] values are converted into the 
metal content Z by means of a relation that allows for 
a-enhancement at decreasing [Fe/H], as suggested by 
Fuhrmann (1998) data. At any age, [Fe/H] was assumed 
to have a la dispersion of 0.2 dex. 

The metallicity of the halo stars was assumed he Z = 
0.0095, with a dispersion of 1.0 dex. This was based on 
an observed [Fe/H] value of — 1.6± 1.0 (Henry & Worthey 
1999), allowing for an a-enhancement of 0.3 dex. 

3.5. The intrinsic colour-magnitude diagrams and 
luminosity functions 

At this point, having described the default SFR, AMR, 
IMF, and libraries of evolutionary tracks and spectra, it 
is useful to open a parenthesis and illustrate the intrin- 
sic CMDs and LFs that we derive from these ingredients. 
Notice that these intrinsic data are the only stellar ingre- 
dients that enter in the old "empirical 4){M) approach" 
mentioned in SecL^ 

In Fig. Owe show, for both disc and halo, the My vs. 
B — V diagram that follow from our choices, together with 
the intrinsic distributions of My and B — V values. 

Several aspects of this figure are remarkable: 

— It is evident that we have a very complete sampling 
of the possible evolutionary stages of the stellar popu- 
lations required in building a realistic Galactic model. 



Limiting to the bright part of the diagrams, we can 
notice the presence of the young disc main sequence, 
subgiants and giants, red clump and horizontal branch 
stars, that make most of the stars observed in shallow 
photometric surveys. 

— The inclusion of white dwarfs, very low-mass stars and 
brown dwarfs allows us to simulate all faint stellar ob- 
jects, including the extremely dim ones - as dim as 
My ~ 44. These faint objects are expected to appear 
in significant numbers in very deep surveys. 

— One should keep in mind that we know, for each sim- 
ulated star, many other physical parameters like the 
complete photometry, age, metallicity, surface gravity, 
etc., so that the same kind of plot could be constructed 
in a multitude of other ways and dissected by grouping 
stars in different parameter bins. 

In summary, the present plots are an evindence of the 
significant advantages of an evolutionary synthesis tool, 
over the empirical approaches which were used in the past 
to construct similar data for Galaxy models. 

Let us also briefly comment on the general aspect of the 
CMDs shown in Fig.|31 In the case of the disc population 
shown in the left panel, it is evident that the low-mass 
and brown dwarf models, taken from different sources, 
combine in a continuous and well-behaved way with the 
sequence drawn by more massive stellar models. There is 
just one abrupt change in the width of this main sequence, 
occurring at My ~ 19, that is caused by the fact that be- 
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Fig. 4. The same as Fig.Obut limited to the disc intrinsic 
Mk vs. J—K diagram. 



low this limit we rely on solar-metallicity models, whereas 
above it the metallicity dispersion is fully represented in 
the models. Exactly the same problem is present for the 
halo population shown in the right panel. Anyway, this 
seems a minor problem because - as will be shown in the 
following sections - these stars, although always present, 
do not make the main features of CMDs observed up to 
now. Also, this seems an acceptable price to pay for having 
an extremely complete intrinsic CMD. 

Finally, we recall that star counts of dwarfs below 
My ~ 6 are affected just by the particular choice of IMF, 
whereas above this limit also the SFR and AMR play a 
major role. 

For the sake of illustration. Fig. 0| shows the intrinsic 
Mk vs. J—K diagram. This looks very different from the 
former i?y-diagram, and is particularly useful for the dis- 
cussion of 2MASS data (see Sect. 15.31 and Marigo et al. 
2003). 



3.6. Geometry of Galaxy components 

Five are the Galaxy components presently defined in 
TRILEGAL: the thin and thick discs, the halo, the bulge, 
and the disc extinction layer. There is also the possibility 
of simulating additional objects of known distance. 



The thin disc: Its density is assumed to decrease exponen- 
tially with the galactocentric radius projected onto the 
plane of the disc, R, 



Pd = Cd exp(-i?//i_R) /(z) 



(6) 



The vertical distribution f(z) is either an exponential, 
exp{—\z\/hz), or a squared hyperbolic secant fimction, 
sec]i^{2 z/hz). Importantly, the vertical scale heigth hz is 
assumed to increase with the stellar age t according to the 
formula suggested by Rana & Basu (1992): 



hd{t) = zo {l+t/tof 



(7) 



where zq, Iq and a are adjustable parameters. This means 
that stars are formed very close to the Galaxy Plane, with 
a scale heigth zq, and later disperse vertically. For practi- 
cal reasons, the code has to deal with a limited number of 
age intervals and scale heigths. Thus, the total age inter- 
val is divided into a number of subintervals - at least 
10, but typically 100 - and the stellar densities computed 
separately for them. 

The normalization constant Cd is set so as to produce 
a given "total surface density of thin disc stars ever formed 
in the Solar Neighbourhood" , 



Sd(0) 



+ 00 



Pd(r)i?=_R0 dz 



where V'd {t) is the SFR per unit disc area in the so-called 
Solar Cylinder^, and Iq is the Galaxy age. So, the thin disc 
geometry is completely defined by the parameters I]d(©), 
^R: Zq, to and a. Additionally, the parameter may be 
adjusted to provide better accuracy at the cost of larger 
CPU times. In present applications we use A'd as large as 
100, which in fact represents an excellent accuracy. 

The thick disc: Similarly to the thin disc, it is described ei- 
ther by a double exponential or by an exponential times a 
sech^ function. However, the scale height is assumed to be 
independent of age. This because the thick disc is always 
described by predominantly old {t > 10 Gyr) populations. 

So, just the parameters l]td(0), ft-_R,,td, and hzM, de- 
fined similarly to the thin disc ones, would suffice to define 
the thick disc. Needless to say, this galactic component 
can be incorporated in the formula for the thin disc, by 
assuming suitable scale heights at large ages in Eq. (0 . 

The halo: Its density ph is given by a de Vaucoulers 
(1959) r^^^ law, deprojected according to Young (1976). 
Alternatively, an oblate spheroid (Gilmorc 1984) can be 
assumed. The halo parameters are: the radial scale rh, the 
oblateness 6h, and the local "total density of halo stars 
ever formed in the Solar Neighbourhood" , 

Jo 



The Solar Cylinder refers to a cylinder perpendicular to the 
Galaxy Plane, of small diameter and infinite height, centered 
on the Sun. 
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The bulge: Although it is included in the code as a triaxial 
truncated spheroid, it has not yet been calibrated. For 
this reason no bulge component is included in any of the 
calculations presented in this paper. 

The disc extinction layer: It is assumed to have an 
exponentially-decreasing density in the vertical direction, 
with a scale-height h'^'^^^ (Parenago 1945, cf. Mendez & 
van Altena 1998). The increase of Ay with distance is 
proportional to this density. This geometric distribution 
is normalized in two possible ways: either adopting a lo- 
cal absorption density Ay, of about 0.75 mag/kpc (Lynga 
1982), or adopting a total absorption at infinity Ay as 
given by Schlegel et al. (1998) maps. We set the second 
option to be the default one. We adopt ft.^""^' = 110 pc as 
Lynga (1982). 

Extinction is always specified in terms of Ay {V in this 
case stands for the effective wavelength of Johnson's V- 
band, 5550 A). The Ay values found for individual stars 
are later converted to those in the several pass-bands, ^a, 
using the Ax/ Ay ratios previously tabulated fSect. lT^ . 

Additional objects: They can be inserted at known dis- 
tance, and optionally assuming a known foreground ab- 
sorption. This alternative is useful for including objects 
such as star clusters or nearby satellite galaxies in the 
simulations. Of course, specifying the SFR and AMR is 
also necessary, which for star clusters is limited to a single 
age and metallicity value. The object total mass is again 
expressed in units of total numbers of stars ever born in 
that area of the sky. 

This option was recently used by Marigo et al. (2003) 
to simulate the 2MASS data towards the Large Magellanic 
Cloud (HMC), and by Carraro et al. (2002) to simulate the 
field of the open cluster NGC 2158. Open cluster simula- 
tions in EIS fields will also be the subject of an upcoming 
paper (Hatziminaoglou et al., in preparation). 

3. 7. Additional parameters 

We also have "pointing parameters" , that specify the re- 
gion of the Galaxy sampled during a simulation. Two main 
modes are presently allowed: either (1) simulations of a 
projected (conic) region of the sky, that requires the spec- 
ification of the central Galactic coordinates {£, b) and total 
sky area, or (2) a volume-limited sample centered on the 
Sun and complete up to a specified maximum distance. 

In both of the above cases, by specifying a given lim- 
iting magnitude in any of the available filters, we avoid 
generating too many faint stars in the output. 

Another input parameter is the resolution in magni- 
tudes. Am. It represents the largest sub-step for the nu- 
merical integration of Eq. Any detail of the Galaxy 
geometry that is caused by a depth structure (in distance 
modulus) smaller than this resolution, or by a LF struc- 
ture finer than it, will be lost. A resolution of 0.1 mag is 
adequate for the purpose of this paper. 



The Sun's position with respect to the Galactic Plane 
is specified by the galactocentric radius and height on disc. 

4. The initial calibration 

The initial calibration is described in Paper I. The most 
important points, relevant for the present paper, are re- 
peated here. That paper describes, amongst other things, 
the first application of the TRILEGAL code. The initial 
calibration is derived from the six fields at high galactic 
latitude covered by the "Deep Multicolor Survey" (DMS; 
Osmer et al. 1998, and references therein), and EIS data 
for the South Galactic Pole (SGP; Prandoni et al. 1999). 
Then, the code, with the parameters fixed, was applied to 
the EIS data in the Chandra Deep Field South (CDFS; 
Arnouts et al. 2001, Vandame et al. 2001). 

4.1. Paper I calibration 

The IMF, SFR, and AMR for the disc and halo were those 
already specified in Sects. rOI and as being the default 
ones. 

The disc component was described by a double- 
exponential in scale height and Galactocentric distance. 
The model did not have separate components represent- 
ing the thin or thick disc. Instead the scale height for disc 
stars was a function of age, and was parametrized as in 
Eq. Q. The parameter values in this equation were not 
the same as in Rana & Basu (1992) - namely zq = 95 pc, 
to = 0.5 Gyr and a = (2/3), since this does not fit very 
well the derived scale height of "thick" , "old" , "interme- 
diate" and "young" disc components as derived by Ng 
et al. (1997). Their results are described by zq — 95 pc, 
Iq — 4.4 Gyr and a — 1.66, which was adopted in Paper I. 

Since none of the six DMS fields, nor the CDFS and 
SGP, contains a bulge component, this population was not 
included. 

The Sun was assumed to be 15 pc above the Galactic 
Plane (Cohen 1995; Ng et al. 1997; Binney et al. 1997) 
and the distance of the Sun to the Galactic Centre was 
assumed to be 8.5 kpc. 

With these input ingredients fixed, the halo oblateness 
(and local halo number density) was derived by fitting the 
number of halo stars, defined by (in the Johnson-Cousins 
system) < B-V < 0.7 in the range 20 < S < 22 and 
< V-I < 0.8 in the range 18 < / < 20, in these seven 
fields and was found to be q = 0.65 ± 0.05. This value 
was smaller than the value of 0.8 ± 0.05 quoted by Reid & 
Majewski (1993), but Robin et al. (2000) could not exclude 
a spheroid with a flattening as small as q = 0.6 and Chen 
et al. (2001) derived q = 0.55 ± 0.06. 

The disc radial scale length (and local disc number 
density) was derived by fitting the number of disc stars, 
defined by 1.3 < B-V < 2.0 in the range 20 < B < 22 
and 1.8 < V-I < 4.0 in the range 18 < / < 20, and was 
found to be hn — 2800 ± 250 pc. This was in agreement 
with the lower limit of 2.5 kpc (Bahcall & Soneira 1984) 
and the work of Zheng et al. (2001) on M-dwarfs who 
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Fig. 5. Distribution of stellar distance moduli in the sim- 
ulation corresponding to our initial calibration, in a conic 
bean towards the NGP. This is shown for 11 disc com- 
ponents of increasing age - at steps of 1 Gyr, as labelled 
- and for the halo. The top panel shows all the curves, 
whereas the bottom one expands the vertical axis in order 
to detail the profiles for ages younger than 5 Gyr. It can be 
noticed that younger disc components are found at lower 
mean distances (peaking from say = 8 to 12 as the age 
increase) , whereas the halo stars are found with a nearly- 
Gaussian distribution of /xq which peaks at about 15. The 
characteristic shapes of these distributions result from two 
competing trends in Eq. the quadratic increase of the 
geometrical factor , and the almost-exponential decrease 
of the density p with distance to the Galactic Plane. 



disc, we separate the profiles coming from different ages 
(i.e. different scale heights) spaced by 1 Gyr. As can be 
noticed, a simulation of the NGP - if not constrained by 
any limiting magnitude - would contain increasing num- 
bers of disc stars as we go to older ages, and at increasing 
mean distances (from /iq ~ 6 to 12 as the age goes from 
very young to 11 Gyr). For each age considered, the disc 
distribution of /iq looks like an asymmetric curve with a 
slow increase followed by a faster decay. The fiQ distribu- 
tion of halo stars, instead, looks like a single Gaussian of 
mean /ig ~ 15. 

As a rule, we can conclude that halo stars dominate 
star counts at very large distance moduli (/^o ^ 13), 
whereas intermediate-age to old disc stars would dominate 
counts at most "intermediate distances" {po from say 9 to 
13). Only at very short distance moduli - /iq ^ 9i i-S- 
in the immediate Solar Neighbourhood - can the young 
disc stars make a sizeable contribution to the star counts. 
Alternatively, one has to look at lower galactic latitudes 
to see a higher contribution from the young disc. 

Of course, the situation gets more complex as we con- 
sider the limiting magnitudes that are present in any sur- 
vey, and that favour the detection of the few closest stars 
in spite of the many distant ones. Anyway, the present 
Fig. Elshows the type of stars which make the major con- 
tribution depending on the depth of a given survey. This 
information is relevant for the discussion presented in the 
next section. 



5. Recalibration and fine tuning 

Since Paper I, we have improved many aspects of 
TRILEGAL, and checked the model predictions with ad- 
ditional datasets. This has forced some changes in the 
model calibration, as will be detailed in the subsections 
below. 



derived Hr = 2750 ± 160 pc and Ojha (2001) who derive 
hn = 2800 ± 300 pc for the thin disc. 

The model with these parameters was then used to 
estimate the stellar counts in the CDFS field, yielding a 
fairly good fit of the UBVRIJK number counts, CMDs, 
and colour distributions. The model with these parame- 
ters was also used by Marigo et al. (2003) to successfully 
predict the foreground population towards the LMC in 
JHK. 

4.2. Distribution of distance moduli 

It is very instructive to look at the characteristic dis- 
tributions of distance moduli, fiQ = 51ogr — 5, for this 
model calibration. Were all stars in a given field - even the 
dimmest ones - possible to be observed, such distribution 
would be proportional to the integral of the quantity p 
(see Eq.nj over small bins of distance modulus. We show 
these quantities as evaluated for the line of sight of the 
North Galactic Pole (NGP) in Fig. In the case of the 



5.1. Summary of recent changes 

First of all, we opted for a better temporal accuracy as 
given by = 100 (see Sect. OSll instead of the = 10 
adopted in Paper I. = 100 implies a virtually continu- 
ous change of thin disc geometry with age, which in itself 
represents a novelty in star count models. 

Second, we have adopted a more realistic metallicity 
distribution for halo stars, namely the one measured by 
Ryan & Norris (1991). The Sun's height above the disc 
was corrected to 24.2 pc (cf. Mai'z-Apellaniz 2001). 

So far, these changes have just a minor impact in our 
final number counts. Major revisions instead resulted from 
our attempts to reproduce 2MASS and Hipparcos data, 
not considered in Paper I. Without entering in a detailed 
description of all attempts carried out to fit the models 
to the data, here we describe the main arguments used in 
establishing our final calibration. They are: 

1) The Paper I calibration, although in excellent agree- 
ment with deep star counts and with the constraints from 
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Fig. 6. Number counts as a function of magnitude. The histogram with error bars represents the data from the 7 and 
5-passband CDFS steUar catalogues (Paper I) with poissonian errors. The dashed and continuous gray vertical lines 
indicate the magnitude limits for efficient morphological classification, MAG_STAR_LIM, and the 90 percent completeness 
limit, respectively. For the J and K plots, at brighter magnitudes (J < 16.5 and K < 15, respectively) we add 
the histogram corresponding to 2MASS JKs data. The smooth lines refer to the model results for a region of much 
larger area and hence better statistics, computed with the same binning (0.5 mag) as the data: they show separately 
the contribution of the halo (dot-dashed line), of the disc (dashed), and the total result (continuous line). At the 
top of the plot, we present the number ratio between observed and modelled sources as counted within the limits 
indicated by short gray vertical marks (the fainter limit coincides with MAG_STAR_LIM). Inside these magnitude limits, 
the comparison between the continuous line and histogram indicates the goodness of the model in reproducing the 
observed star counts. 

CDFS 220.0 -53.9, -0 4<U-B<1.2 CDFS 220.0 -53.9, -0.4<U-B<1.2 




Fig. 7. The same as Fig.El but now limiting data and models to the "blue subsample" with —0.4 < U — B < 1.2, which 
is dominated by the disc at brighter magnitudes and by the halo at the faintest. To help in the comparison, the scales 
are kept the same as in Fig. The comparison between data (histogram with error bars) and models (continuous 
line), inside the magnitude interval delimited by the vertical gray marks, probes whether the models correctly predict 
the number ratio between halo and disc stars (shown separately by dot-dashed and dashed lines, respectively). 
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Hipparcos, presented systematically too few stars at in- 
termediate magnitudes, of say K ^ 14, as compared to 
2MASS data. The deficit was small for high latitudes, 
reaching factors of about 2 at 6 ^ 10 deg. This calls for an 
decrease in the slope of the LFs. 

2) Star counts at faint magnitudes are determined 
essentially by the halo plus the oldest disc components 
(those with highest scale heights, see Fig. EJ. A way to 
decrease the LP slope is to reduce the contribution of the 
oldest disc components, and then compensate the reduc- 
tion in total number counts by increasing the total disc 
SFR. We perform such an operation in the following way: 
the age scale of the thin disc SFR is changed by a constant 
factor of 0.8 (or —0.1 dex in logi), so that the oldest disc 
component is found at an age of 9 Gyr, and with scale 
height of hz — 603 pc. 

3) This change in the age (and height) scales of the disc 
is quite satisfactory in predicting star counts in CDFS, 
DMS, and 2MASS, but now fails by predicting twice as 
many stars as sampled by Hipparcos. This local sample 
is dominated by relatively younger stars, of ages ^ 2.5 
Gyr (as indicated by Fig.EJ and ^ 150 pc. In order to 
eliminate this discrepancy, we simply change the vertical 
distribution of thin disc stars, from the usual exponential 
to a squared hyperbolic secant law with twice its scale 
height, i.e. 

exp{-\z\/h,) 0.25 sech^ {0.5 z/h,) (8) 

In this way, star counts at deep and intermediate mag- 
nitudes do not change by much, while the local one is 
reduced to the level of the observed values. 

The assumption of a sech^ function instead of an expo- 
nential one is far from being arbitrary. Exponential laws 
are commonly assumed in star count models mainly be- 
cause this is inferred from the surface brightness profiles 
in the discs of edge-on galaxies. But this observation does 
not refer to the inner part of the discs, which are nor- 
mally obscured by dust lanes and are the site of ongoing 
star-formation. From the kinematic point of view, theory 
predicts that the density profile of isothermal discs do fol- 
low a sech^ law - then approaching an exponential law 
at high z, in accordance with observations. Of course, the 
disc of our Galaxy is not isothermal - as demonstrated by 
the observed increase of scale heights with the stellar ages 
-, but this latter theoretical aspect makes the sech^ law 
to be a better approximation than a simple exponential 
law. 

Finally, in order to improve upon our final results for 
2MASS and Hipparcos, he have modified the SFR of the 
thin disc: we have assumed that between 1 and 4 Gyr the 
SFR has been 1.5 times larger than at other ages. This 
causes a moderate impact in the distribution of stars in 
the My vs. B — V diagram, that can only be explored by 
means of Hipparcos data. In fact, such a change in the 
disc SFR has a negligible impact on deep fields and just a 
minor impact on 2MASS counts. 

The normalization constants we derived imply 



— a local surface density of ever-formed disc stars of 
Sd(0) — 59 M0pc~^, a number that compares well 
with the present dynamical surface density of matter 
in the Solar Cylinder (56 ± 6 Mq pc~^ cf. Holmbcrg & 
Flynn 2004). We find reassuring that these two quan- 
tities, being very different in their nature, have the 
same order of magnitude. This would indicate an effi- 
cient conversion of baryonic matter into stars over the 
history of our disc. 

— a local volume density of halo stars ever-formed 
17h(0) = 1.5 X 10-4 M0PC-3. 

Let us now present the results for this calibration, pro- 
viding more details on the way the different data samples 
have been selected and modelled. 

5.2. Simulating deep fields: CDFS, DMS, and SGP 

A main characteristic of deep photometric surveys is the 
rich presence of galaxies, a significant fraction of which ap- 
pears as point sources and cannot be easily distinguished 
from real stars. Thus, object classification by means of 
morphological and photometric criteria is of central im- 
portance in these fields. The reader is referred to Paper I 
and Hatziminaoglou et al. (2002) for a discussion of these 
aspects. 

In this paper, we deal with 3 deep catalogues which, 
as far as possible, have been cleaned from contamination 
by galaxies. 

5.2.1. CDFS 

The first one is the CDFS stellar catalogue (Paper I), 
which points towards a relatively clean area centered 
in {£ = 220.0,6 = -53.9). According to Schlegel et 
al. (1998) maps, the reddening for background sources 
amounts to Eb-v = 0.0148 or Ay = 0.0458. We make 
use of UBVRI data from the 5-passband catalogue cov- 
ering 0.263 deg^, and JK from the 7-passband catalogue 
covering 0.0927 deg^. The data has been cleaned from non- 
stellar objects according to the criteria and methods thor- 
oughly discussed in Paper I. For this field, we simulate the 
galactic population for a 2.63 deg^ region. Fig.^lshows the 
results in units of number of stars per unit deg^ and 0.5- 
mag intervals. As can be noticed, the agreement between 
modelled and observed counts is good, their ratio being 
comprised between 0.87 and 1.08 for all 7 pass-bands con- 
sidered. 

Particularly interesting is the comparison between sim- 
ulated and observed number counts for stars in the inter- 
val —0.4 < U — B < 1.2. This subsample, according to the 
models, is completely dominated by disc stars at brighter 
magnitudes {V ^ 18), and then by halo stars at fainter 
magnitudes {V ^ 20); this is illustrated by the compar- 
isons between observed and simulated CMDs shown in 
Paper I (see in particular the online version of their fig- 
ure 5). Hence, by fitting number counts in this particu- 
lar subsample we can be sure to be correctly modelling 
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the relative proportion between halo and disc densities. 
Moreover, a "blue subsamplc" defined in this way is com- 
posed mostly by stars in well-modelled evolutionary stages 
- i.e. main sequence stars of moderately high Tcff" and ex- 
cludes the hot white dwarfs and the reddest very-low mass 
stars, for which the reliability of present-day evolutionary 
and spectral models could still be questioned. 

Such a comparison is presented in Fig.d As can be no- 
ticed, the blue subsample presents good evidence that the 
ratio between halo and disc densities is well represented 
for this line-of-sight. The most significant discrepancy that 
appears at this point is a moderate excess of halo stars at 
B ~ 23. As demonstrated in Paper I, this is a magnitude 
interval in which number counts become sensitive to the 
very low-mass IMF. In fact, we have verified that a bet- 
ter agreement with observations turns out if we artificially 
eliminate halo stars with Mj ^ 0.2 Mq from our models. 
However, we prefer not to draw any strong conclusion from 
this test, since ~ 23 is also close to both the limit for 
efficient morphological classification, MAG_STAR_LIM, and 
the 90-percent completeness limit of CDFS data. 

5.2.2. DMS 

The Deep Multicolor Survey by Hall et al. (1996) and 
Osmer et al. (1998) provides deep UBVR'm and 786 
data for 6 different fields of |6| ^ 30deg. From their cata- 
logue, we eliminated the objects that are likely not to cor- 
respond to stars, i.e. those marked as "galaxy" , "noise" , 
"diffuse object" , or "long object" in any of the DMS pass- 
bands. 

Figure IHl shows the UBVR results as compared to our 
models. As shown in the figure, if we limit the analysis to 
magnitudes between the DMS "upper limiting" and the 
"threshold" one - an interval almost free from saturation, 
incompleteness and contamination effects - there is a good 
overall agreement in the number counts between the data 
and our model. With the exception of DMS field 21, the 
largest discrepancies reach about 10 percent in number 
counts. 

Again, the comparison between simulated and ob- 
served number counts for stars in a "blue subsample" - 
now defined in the interval —0.4 < U~B < 0.8* - are 
important for verifying whether we have the right propor- 
tions between disc and halo stars. A careful examination 
of Fig. [^reveals that this is indeed the case, but for DMS 
field 21, where, apparently, simulations contain too few 
disc stars. Field 21 is the innermost of DMS fields, pointing 
to a Galactic region for which actually our model calibra- 
tion indeed presents the largest problems (see Sect. I5.3|l . 



* The Johnson U — B colour of DMS separates less the red 
and blue sequences typical of field CMDs, than the U — B 
colour of CDFS - this latter being based on ESO WFI fil- 
ters. This is the reason why we have used a smaller interval 
of U — B to define the blue subsample of DMS data, than for 
CDFS. 



5.2.3. SGP 

The South Galactic Pole (SGP) as observed by EIS 
(Prandoni et al. 1999; Zaggia et al. 1999) presents BVI 
data for an effective area of 1.21 deg^ centered at £ ~ 
306.7 deg, 6 = — 87.9deg. Fig. [TUI shows the comparison 
between the data and our models. This time, the shape 
of the observed counts are, in general, fairly well repro- 
duced, except that the amplitude of our model is twice 
as large as observed. We did not find any obvious way of 
reducing this discrepancy, without spoiling the excellent 
agreement we find for the other fields considered in this 
paper - including 2MASS data for the SGP itself. Also, 
the lack of U -band data prevent us from analysing a blue 
subsample essential in order to investigate the source of 
the discrepancies. 

5.3. Simulating 2MASS data 

As seen above and also in Paper I, the deep simulations 
of DMS, SGP and CDFS data are ideal for probing the 
relative proportions between halo and disc components, as 
well as the shape of the halo and the IMF of the disc. For 
better probing the disc and its details (spiral arms, dust 
lanes, warps, etc), shallower surveys covering larger areas 
are better suited. However, if one wants to avoid compli- 
cations caused by dust, we should consider only counts at 
the infrared. In this context, 2MASS constitutes an invalu- 
able dataset: it covers the all sky in JHKs for magnitudes 
as faint as J ~ 17. 

From the 2MASS All-Sky Data Release Point Source 
Catalog we have selected the sub-sample obeying the so- 
called "2MASS level 1 science requirements" (see the 
User's Guide in Cutri et al. 2003). In practice, these crite- 
ria refers to stellar sources falling outside of tile overlap re- 
gions, and of high photometric quality (namely S/N > 10 
and a < 0.11 mag, band- by-band) . For most of the sky - 
excluding the most crowded low-latitude and bulge fields - 
this subsample of 2MASS is essentially complete for mag- 
nitudes brighter than about 15. 

In the panels of Fig. 1111 we show the complete results 
for two 2MASS fields, one of high latitude (the NGP at 
e = deg, 5 = 90 deg) and one of low {£ = 180 deg, b = 
lOdeg). Plots for symmetric fields - namely the SGP and 
the {£ ~ 180 deg, 6 = —10 deg) ones - look very much 
the same and present similar number counts. The counts 
in these particular fields are very well reproduced by the 
model, over a wide range of magnitudes in all the 3 pass- 
bands of 2MASS. The reader can also notice that the J, H 
and Kg diagrams contain essentially the same information, 
so that examining all of them may be redundant. For this 
reason and for the sake of conciseness, in the next figure 
we prefer to present just results for the H band. 

Fig. El presents the i?-band results for a series of 
2MASS fields disposed along a great circle in the sky - the 
one ai £ — 180 deg, encompassing the Galactic centre, 
anticentre and polar regions. We show only the northern 




Fig. 8. The same as Fig. El but for the several fields of the DMS (Osmer at al. 1998), whose {i,b) are given at the 
top of each panel. The data (histogram with error bars) excludes the objects that are likely not to correspond to 



Girardi et al.: Star counts 



15 






Fig. 9. The same as Fig.|Hl but now limiting data and models to the "blue subsample" with —0.4 < U — B < 0.8. 
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Fig. 10. The same as Fig. El but for the ElS-deep SGP 
data (Prandoni et al. 1999). The hmits of rehabihty of 
the data were located, somewhat arbitrarily this time, at 
magnitudes 16 and 21.5 for all filters. 



Galaxy fields, since results for southern fields are essen- 
tially the same. 

The results are certainly encouraging. We are able to 
predict the correct number counts, with errors smaller 
than ~ 30 percent, for all fields located at least 10 de- 
grees above the Galactic Plane, except for inner Galactic 
regions where the lack of a Bulge component - presently 
not included in the model - becomes noticeable. Moreover, 
it is not to be excluded that the present description of the 
Galactic halo fails for small galactocentric distances, hence 
contributing to the discrepancies we find at inner galactic 
regions. 

Remarkably, we find that for all 2MASS fields away 
from the Bulge we have analysed, disc stars make the bulk 
of the number counts. The halo contribution is almost 
negligible. However, we recall that halo stars contribute 
to make a particular feature in wide-area J vs. J — Ks 
CMDs, namely the central of three vertical fingers dis- 
cussed by Marigo et al. (2003), that are clearly present in 
2MASS CMDs towards the LMC. This indicates that the 
halo component has its importance in analyses of 2MASS 
photometry. 

In fields too close to the Galactic Plane, instead, our 
model predictions fail, as demonstrated by the first panel 
of Fig. El regarding the direction of the Galactic anti- 
centre. One of the reasons for this failure is surely the 
too simplistic modelling of the dust distribution along the 
Galactic Plane. Improving this aspect of the model, how- 
ever, is beyond the scope of the present paper. 



5.4. Simulating Hipparcos data 

Examining the stellar counts in the immediate Solar 
Neighbourhood is an obvious test for any Galaxy model. 
In fact, present models usually check whether their results 
are consistent with the observed "local stellar density" , or 
with some similar parameter derived from Hipparcos data 
(see e.g. Robin et al. 2003). Of course, by using just a 
single density parameter as a constraint, we ignore the 
wealth of photometric information that is present in the 
data for nearby stars, that could tell us much about the 
distribution of stellar parameters in the disc. In order to 
start exploring this information, in the following we try to 
predict the counts of a local sample, using the Hipparcos 
dataset. 

First of all, however, we should remind that our model 
in Paper I has been effectively calibrated on deep data 
{V ^ 15). For such deep surveys the effective distribu- 
tions of distance moduli should approach the ones shown 
in Fig. i.e. the samples are dominated by relatively far 
objects, at distance moduli ranging from /xq ~ 7 to 14 for 
the disc populations, and from /io ~ 12 to 18 for the halo. 
It is then obvious that our previous model calibration has 
almost nothing to do with the very local sample of stars, 
i.e. the one observed at distances lower or comparable to 
100 pc, and at bright magnitudes such as y ^ 8, that we 
will define below. Deep and local samples could even be 
considered, in terms of their stellar populations, as com- 
pletely independent data samples. 

The Hipparcos and Tycho catalogues (ESA 1997) have 
provided parallaxes with ~ 10 milliarcsec (mas) accu- 
racy for several thousands of stars, together with accu- 
rate BV photometry. The Hipparcos input catalogue was 
constructed in such a way that there are no clearcut cri- 
teria for defining volume-limited samples out of its data. 
This problem has been recognized by a number of previ- 
ous authors. Hernandez et al. (2000) and Bertelli & Nasi 
(2001), for instance, find it to be extremely difficult to 
define volume- limited samples containing enough stars for 
studying the SFR in the solar vicinity up to the oldest 
possible ages (see also Schroder & Pagel 2003). 

However, our aim in this paper is different from pre- 
vious works. We consider a subsample of the Hipparcos 
catalogue to be good provided it is complete and could be 
simulated. Differently from the above-mentioned papers, 
we do not need to limit our simulations to stars being all 
contained in the same volume. Our sample can be selected 
by using a few simple criteria, based on the following re- 
alizations: 

— the Hipparcos catalogue contains all stars brighter 
than Vvnn ^ 7; 

— from them, all parallaxes tt higher than 7rii,„ ~ 10 mas 
have been measured; 

— due to the presence of parallax errors, the sample de- 
fined by TTiini comprehends a volume of radius riim 
somewhat larger than l/7riim — 100 pc; 

— among these stars, the bulk of binaries has been iden- 
tified. 
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Fig. 11. The same as Fig. but for two sample fields of 2MASS data: {£ = Odeg,6 = +90deg) (the NGP), and 
{£ = 180 deg, b — +10deg). Lines have the same meaning as in previous figures. In both cases, the disc component is 
responsible for the bulk of the number counts in the model. 



Therefore, we select from the Hipparcos catalogue all stars 
with V < Viiai and tt > nnai, and not classified as binaries 
of any kind. Vnm and TTiim are kept as parameters. Initially, 
we have conservatively adopted Vum = 7 and Triim = 10 
mas, which corresponds to > 100 pc. 

We then simulated this local sample using the 
TRILEGAL code. To do so, we have generated synthetic 
samples up to a distance rii,„ = 200 pc. This is large 
enough to include all stars that, due to parallax errors, 
will later be scattered to apparent distances closer than 
100 pc. 

The simulated true physical distances tq are first con- 
verted in the true parallax ttq = l/?'o- The simulated 
parallax errors Sir (described in the Appendix) are then 
added to ttq so as to generate the "observed" parallaxes 
and distances, tt and r = 1/tt. The "observed" abso- 
lute magnitude is then derived by the usual formula, 
My = V — Slogr + 5. Extinction has been ignored, since 
its effect inside a radius of 200 pc is neghgible. 

The results of this exercise can be seen in Fig. ^1 for 
the conservative choice of V^im = 7 and TTiim = 10 mas. 
In the left-hand panel, we compare the simulated and ob- 
served My VS. B — V diagram. The agreement between 
simulated and observed samples is striking. It can be no- 
ticed that models describe very well both the location and 
width of the main sequence all along from My ~ 9 to —4. 
Particularly good is the description of the left boundary of 
the MS, very well marked because it is produced by stars 
in their phase of slowest evolution, when they depart from 
their ZAMS to the right in the HR diagram. Regarding the 
MS width, we know from stellar models that it is affected 
essentially by two factors: the assumed metallicity disper- 



sion (or equivalently the AMR in our models) and the 
efficiency of convective core overshooting for M > \ Mq 
stars. The good agreement between models and simula- 
tions seems to indicate that both these ingredients are 
well described in our models. Of course, before consider- 
ing the implications that this result may have for the disc 
AMR and for the theory of stellar evolution more quanti- 
tative comparisons would have to be made. However, such 
a discussion is beyond the scope of the present paper. 

Regarding the evolved stars (subgiants and red giants), 
the agreement is also very good: we can notice the right 
width of the subgiant and lowest part of RGB; the clump- 
ing of core-He burning stars in the right location; the bi- 
furcation of the red giants above the clump into two loose 
broad sequences: a vertical one made of intermediate-mass 
core-He burners and the inclined one, going to the red, 
made of genuine first-ascent RGB and early-AGB stars. 

As can be seen, the simulation predicts about the same 
star counts as observed: the total number of objects in 
both panels is 4085 and 4182 for the data and models, 
respectively. 

The right panels show the corresponding distributions 
of apparent distance r and derived absolute magnitude 
My (continuous lines). The dashed hues refer to the sub- 
sample of subgiants and giants, defined by the stars with 
Mv > 6.82 X (B-V) - 2. 

There, although a first look indicates a good overall 
level of agreement between models and data, some dis- 
crepancies become apparent. In the histogram of r, the 
most noticeable one seems to be a modest excess of sim- 
ulated stars at r ~ 80 pc, that amounts to about 20 per- 
cent. Since the volume sampled in the simulation is very 
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Fig. 12. The same as Fig. ^2 but limited to the iJ-band and for a series of fields disposed along the £ = 0, 180 deg 
great circle in the sky. Just the northern part of the circle is presented, since the results are similar for the southern 
fields. For the fields pointing towards the outskirts of the galactic bulge, {i = 0,b — +20) and {i = 0,b = +10), the 
bulge population is clearly seen as an increase in stellar counts at H > 8, caused by the bulge RGB, which is not 
accounted for by our disc+halo model. 



small, we consider that such a discrepancy is unlikely to 
be derived from inhomogeneities in the local distribution 
of stars; more likely, a better simulation of Hipparcos par- 
allax errors could improve the models in this particular 
range of distances (parallaxes). There is also a modest 
deficit of simulated stars in the smallest distances (from 
15 to 35 pc), amounting again to about 20 percent, which 
however does not appear among the subgiants and giants. 
This again might indicate a problem in the simulation of 



parallax errors for the faintest stars. On the other hand, 
the spike of observed star counts at r = 45 pc is to be 
assigned to the Hyades cluster, which was not considered 
in our simulation. 

The histogram of number counts against My indicates, 
once again, some modest discrepancies, that are however 
statistically very significant. The most important one con- 
sists on an excess of simulated bright stars accompanied 
by a deficit of the faintest ones. We have performed a KS- 
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Fig. 13. The several panels show our Hipparcos simulation (grey) versus the real data (dark). On the left panels, 

we have the My vs. B — V diagram of both samples (observed and simulated) limited to a parallax of tt < 10 mas 
(apparent distance r = I/tt < 100 pc) and an apparent magnitude of < 7. The agreement between simulated and 
observed samples is remarkable for most regions of the diagram. The lines in the right panels show the corresponding 
distributions of apparent distance r (top) and derived absolute magnitude (bottom), for both the total samples 
(continuous lines) and the subsample of subgiants and giants (dotted lines). For the data histograms, error bars 
denote the standard error {^/N) of a Poisson distribution, and serve as a guide to the comparison between model 
and data. As can be noticed, the simulation predicts about the same total star counts as observed, with just a 2.5 
percent excess of model stars. Looking at the right panels, however, we notice a significant excess of simulated stars at 
r ~ 80 pc, a deficiency at r ~ 20 pc, and the obvious failure to reproduce the spike at r = 45 pc that corresponds to 
the Hyades cluster. Moreover, in the comparison of number counts as a function of My, there seems to be an excess 
of bright stars and a deficit of faint ones. 



test comparison between the two My distributions, and 
find that the probability of them being drawed from the 
same distribution would be highly increased - from its 
present close-to-zero value up to about 0.3 - if our models 
were shifted by 0.26 mag in My. We think however that 
applying such a shift whould not be justified, and it could 
not be easily translated onto a physical interpretation (i.e. 
shifts in the My distribution can be forced by using mod- 
els with a corrected photometric zero-point, different IMF, 
different SFR, modified prescription for simulating paral- 
lax errors, or a combination of all these effects). 

To conclude, we remark that all the above-mentioned 
discrepancies are, from the statistical point of view, highly 
significant, since they refer to samples containing large 
numbers of stars. They may be indicating points where 
the present models can be improved. They however do not 
invalidate the present model calibration, for a series of rea- 
sons: first of all, many of the discrepancies are suspected 
to result from the imperfect simulation of parallax errors; 
second, when one deals with normal star counts data, fine 
details (i.e. those seen at ^ 0.25 mag resolution) of the 



My distribution become of little relevance since stars are 
dispersed over a large range of distance moduli; third, the 
most relevant comparison regards the total star counts 
that are directly linked to the local density of stars: in 
our case they are very well predicted, to within 5 per- 
cent. Reaching such a result for the star counts is already 
remarkable, whereas reaching a statistically-robust com- 
parison with Hipparcos, in all its details, may be still far 
from reach. In this regard, the present work represents just 
the first attemp. 

Figure presents the same as Fig. El but for a 
slightly deeper sample, of Vum = 8 and Triim = 10 mas. 
The effect of selecting a deeper Viim is that more stars 
with My > 2 are included in the sample, then increas- 
ing the contribution of intermediate-age to old stars (in 
their main sequence, subgiant branch, and initial RGB 
evolution) to the number counts. In this case, the star 
counts are 8055 and 7640 for data and models, respec- 
tively. The discrepancies between model and observations 
seem slightly increased, as expected since we are including 
data for which the completeness starts to become an is- 
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Fig. 14. The same as Fig. but for samples limited to an apparent magnitude of l^i,„ < 8. Notice the increased 
number counts of stars with My > 2. 



sue, and for which parallax errors are slightly larger with 
respect to the previous V < 7 sample. 

6. Discussion and conclusions 

We have presented a new code for simulating the photom- 
etry of Galaxy fields. 

The code has been calibrated by predicting counts 
in a variety of stellar surveys, that comprehend some 
very deep multi-passband catalogues cleaned from galax- 
ies (CDFS, DMS, SGP), the "intermediate-depth" near- 
IR point source catalogue of 2MASS, and the very local 
stellar sample derived from Hipparcos catalogue. 

The results are certainly satisfactory, since we have 
demonstrated that the predicted star counts agree well 
with the observed ones. The typical discrepancies are 
smaller than 30 percent for most of the sky, and inside the 
estimated magnitude limits of reliability of the observed 
star counts. This agreement is remarkable, when we con- 
sider the wide ranges of magnitudes, passbands (from U 
to K) and sky positions that were considered in this work. 

The major discrepancies were found for: (1) Inner 
Galactic fields, located within about 30 degrees from the 
Galactic Center, for which we largely underestimate the 
number counts. Part of this discrepancy can be attributed 
to the lack of a bulge component in the present model, but 
probably a better modelling of the inner disc and halo 
densities is also necessary. (2) Low-latitude fields, with 
|6| ^ lOdeg, which probably require a detailed modelling 
of the distribution of dust and recent star formation along 
the disc. Finally, (3) the SGP as observed by EIS, for 
which we predict twice as many counts in the optical as 



observed. In this case, the origin of the discrepancy could 
not be identified. 

Note that the present model calibration is not yet fully 
optimized and is likely not to be unique, in the sense 
that other choices for the stellar densities and star forma- 
tion histories of Galaxy components might produce sim- 
ilarly good results. In fact, the question arises whether 
TRILEGAL could be adapted to find, in an objective way, 
a maximum-likelihood solution for the functions p, (j), 
etc. - using a database of few high-quality niultiband cat- 
alogues covering several regions of the sky and with a large 
range in depth. The answer is likely yes, but the way is 
certainly not straightforward. Such a work would imply at 
least the following steps: (a) finding a zeroth-order calibra- 
tion producing "acceptable" results, which is actually the 
step performed in this paper; (b) establishing a robust like- 
lihood criterium for comparing the resulting models with 
the stellar data; (c) establishing a numerical algorithm 
to migrate from the zeroth-order to improved maximum- 
likelihood solutions (see for instance Ng et al. 2002); (d) 
exploring the problem of uniqueness of solution by using 
different starting solutions. Therefore, what is presented 
in this paper can be seen as the initial step of a bigger 
project, that we expect to pursue in the future. 

Aside from the good model calibration we have 
reached, the most important advantages of the 
TRILEGAL code can be identified in: (1) the fairly 
complete database of stellar evolutionary tracks already 
implemented; (2) the use of an extended spectral library 
to simulate many different photometric systems, and 
their extinction coefficients, in a self-consistent way, (3) 
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the modular and flexible structure of the code, that 
allows easy changes and additions to both input functions 
(SFR, AMR, IMF, etc) and geometric parameters (the 
density of Galaxy components. Sun's position, pointing 
parameters, etc.). 

With respect to other population synthesis codes com- 
monly used to simulate the photometry of Galaxy fields 
(e.g. Ng et al. 1995, Vallenari et al. 2000, Castellani et al. 
2002, Robin et al. 2003), TRILEGAL shares the advan- 
tage of being intrinsically self-consistent in the relative 
numbers of stars predicted to be in different evolutionary 
phases (including stellar remnants such as white dwarfs). 
In fact, for a given Galaxy geometry, stellar number ratios 
are univoquely determined by the choice of SFR, AMR 
and IMF, and are not tunable parameters. In TRILEGAL, 
this self-consistency of population synthesis codes is kept 
as a very stringent criteria, since there are explicit checks 
for the continuity of all stellar quantities (including core 
mass, envelope mass, and surface chemical composition 
when applicable) in the isochrone-construction routines 
that make part of the code. 

As already mentioned, a main positive characteris- 
tic of TRILEGAL consists in the extreme flexibility in 
the way input libraries (evolutionary tracks, atmospheres) 
and functions (geometry, IMF, SFR, AMR of Galaxy 
components) can be changed, tested, and added to a 
database for future use. Improvements in the stellar evo- 
lutionary tracks, for instance those described in Marigo 
et al. (2003), have been inserted in TRILEGAL in al- 
most no time. We have so far computed test models 
in at least 10 different photometric systems (including 
UBVRIJHK, Washington, HST-based instruments hke 
WFPC2, NICMOS and ACS, the EIS photometric sys- 
tem, 2MASS and SDSS), and we are confident that virtu- 
ally any broadband Vega, AB or ST magnitude system can 
be considered as well (cf. Girardi et al. 2002). Needless to 
say, before the present TRILEGAL calibration has been 
considered as acceptable, we have made wide use of its 
flexibility by testing many different IMFs, AMRs, SFRs, 
and density functions published by different authors. Such 
a flexibility is of fundamental importance for facing the 
huge amount of wide-field photometric data that is be- 
coming available these days, and to take immediate advan- 
tage of the next generation of improved/expanded stellar 
evolutionary and atmospheric models, that are now being 
prepared by different groups around the world. 

Therefore, the TRILEGAL code is ready to use for the 
variety of wide-angle surveys in the optical/infrared that 
will become available in the coming years and will provide 
constraints that will help us to pin down the structure of 
our Galaxy. 
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Appendix A: Simulation of parallax errors in 
Hipparcos 

In this Appendix we describe how the parallax errors in 
the Hipparcos catalogue have been simulated. From the 
about 118000 objects in the Hipparcos catalogue the about 
99000 objects have been selected that fulfill: "goodness-of- 
fit" flag (H30) less than 3, "percentage-of-rejected-data" 
flag (H29) less than 10, "number of components" flag 
(H58) of 1, that have a y-band magnitude and a par- 
allax larger than 0.5 mas. Fig. lA.lj shows the distribution 
of errors cttt for these data. From this dataset the following 
recipe was devised. 

The median parallax error (in mas) is calculated from: 



^mcd 



0.914 ( I 



6.1827 



0.735 



and the minimum parallax error (in mas) from: 

5.6010 



- 0.516 ( —\ 

loy 



0.451 



(A.l) 



(A.2) 



Then a random number is drawn from a Gaussian dis- 
tribution with a mean of 1.0 and a sigma of 0.180. This 
number is multiplied with ct™'''^. If this value is larger than 
cr™'" then this is accepted as the parallax error air. If not, 
(a) new random number(s) is (are) drawn until this is ful- 
filled. 

Once a star of magnitude V and his error cttt are simu- 
lated, the individual measurement error Jtt is drawn from 
the Gaussian of dispersion air. 
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Fig. A.l. Errors in Hipparcos parallaxes. The data is in 
the lower panel, our simulations in the upper one. 



